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ABSTRACT 

This  thesis  gives  a  number  of  approximations  and  bounds 
for  the  renewal  function  in  an  ordinary  renewal  process. 
Each  approximation  and  bound  is  calculated  for  the  uniform, 
gamma  and  hyperexponential  distributions  and  compared  with 
the  renewal  function  for  these  cases.   They  are  also  calcu- 
lated for  the  log-normal  distribution  and  compared  with  re- 
sults of  the  simulation  of  the  renewal  function.   Results 
are  tabulated  and  shown  graphically. 
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I.   INTRODUCTION 

Let  {X, ,X2,...}  be  a  sequence  of  independent  identi- 
cally distributed  non-negative  random  variables.  Let  the 
distribution  function  of  X- ,  i  =  1,  2,...  be  F(x) .   Let 

S   =  X..  +  X0  +    +  X  ,   n>l, 
n    1    2         n*     -   ' 

the  ntn  partial  sum,  and  S  =0.   Define  the  random  vari- 
r  o 

able 

N(t)  =  Sup  {  n  |  Sn  <  t  }, 

which  is  the  number  of  renewals  occurring  before  and  in- 
cluding t.   The  renewal  function,  M(t) ,  is  defined  as 

M(t)  =  E(N(t))  ,    t  >  0.  (1) 

Let  F  ft)  be  the  n-fold  convolution  of  F  with  itself.   The 

n v  J 

relationship  between  M(t)  and  F(t)  is  given  by  the  follow- 
ing, 


ct)  -  f.  y"' 


n=l 

which  leads  to  the  well-known  renewal  equation, 

t 
M(t)  =  F(t)  +  f  M(t  ■  x)  dF(x).  (2) 

o 

The  purpose  of  this  thesis  is  to  investigate  approxima 
tions  and  bounds  for  the  renewal  function,  since  Equations 
(1)  and  (2)  can  only  be  solved  for  special  cases  of  the  dis 
tribution  function  F. 


In  Section  II  formulae  are  collected  for  various  ap- 
proximations and  bounds  that  have  appeared  in  the  litera- 
ture (Barlow  §  Proschan  (1965)  ,  Bartholomew  (1963)  ,  and 
Butterworth  and  Marshall  (1971)).   These  formulae  are 
applied  to  the  uniform,  gamma,  hyperexponential  and  log- 
normal  distributions.   In  Section  III  the  results  of  com- 
putations are  presented  in  both  tabular  and  graphical 
forms  for  each  approximation  applied  to  each  distribution 


II.   APPROXIMATION  FORMULAE 

Throughout  this  thesis  the  following  notation  will  be 
used: 

EtX,]  =  \   , 

A2E[X?] 

k  -  — «-i-  -  1, 


and 


t 

Fe(t)  =  A  J   (1  -  F(u))  du  ,  t  >  0. 
o 

With  this  notation  F   is  an  equilibrium  excess  distribution 

e        ^ 

of  a  renewal  process,  and  from  Smith  (1958),  if  F  is  non- 
lattice  , 


This  equation,  together  with  (2),  plays  an  important  part 

in  the  development  of  approximations  to  M(t) . 

1 .   A  Simple  Lower  Bound 

Butterworth  and  Marshall  (1971)  show  that  for  any 

renewal  process 

n  n 

M(t)>At  +  Ys     Fk™  "  Z  V  Fk-1^  "  V**' 


n  >  0,     (4) 
where  FQ (t)  =  1,  and  when  n=0  both  summations  are  assumed  to 
be  empty.   The  right-hand  side  of  (4)  gives  an  increasing 
sequence  of  lower  bounds  which  converge  monotonically  on 
M(t) .   The  case  n=l  will  be  used  in  this  thesis  and  we  define 

A(t)  =  At  -  Fe(t)  ,  t  >  0.     (5) 


2 .  An  Improved  Lower  Bound 

The  following  general  lower  bound  was  found  by 
Barlow  and  Proschan  (1965)  .   Define 

B(t)  =  p^y  -  1.  (6) 

Then  B(t)  <  M(t) .   Since  B(t)  =  (A(t) /F  (t) ) ,  it  is  clear 
that  B  is  an  improvement  over  A  when  F  (t)  <  1 . 

3.  Upper  Bounds  When  F  is  NBUE 

A  distribution  of  a  non-negative  random  variable  is 
called  NBUE  (new  better  than  used  in  expectation)  if 

00 

f   (l-F(x))        1 

J    l-F(t)    QX  -  X    '  l/J 

t 

for  all  t  where  F(t)  <  1.   Equivalent!/  F  is  NBUE  iff  F  (t) 

>  F(t)  for  all  t.   For  details  concerning  this  and  other 

classes  of  distributions  frequently  used  in  reliability 

theory  see  Marshall  and  Proschan  (1970)  .  Butterworth  and 

Marshall  (1971)  show  that  if  F  is  NBUE, 

n  n 

M(t)  <  At  +  £  Fk(t)  "  Z  VFk-l(t)  '  n  -  °'  (8) 
k=l         k=l 

where  again  F  (t)  =  1,  and  when  n=0,  both  summations  are 
taken  to  be  empty.   The  right-hand-side  of  (6)  gives  a  de- 
creasing sequence  of  upper  bounds  which  converge  monotonic- 
ally  to  M(t)  .  *   The  case  n=2  will  be  used  in  this  thesis. 


i 
If  the  inequality  is  reversed  in  (7) ,  F  is  called  NWUE 
(new  wors ;  than  used  in  expectation).   In  this  case  (8)  gives 
an  improved  sequence  of  lower  bour.ds  over  the  sequence  given 
by  (4). 


We  define 

C(t)  =  At  +  F(t)  -  Fe(t)  (9) 

=  A(t)  +  F(t). 

4 .  Upper  Bound  When  F  is  IFR 

A  distribution  is  said  to  have  increasing  failure 
rate  (is  IFR)  if  Log(l-F(t))  is  concave  in  t.   Under  this  as 
sumption  Barlow  £,   Proschan  (1965)  showed  that  if  we  define 

D(t)  =  **£[*}■  ,  (10) 

then  D(t)  >  M(t) . 

5 .  Approximations  for  M(t) 

Bartholomew  (1963)  derived  the  following  approxima- 
tion for  the  renewal  density  m(t)  =  (dM(t)/dt)  when  F  has  a 
density  f; 


m(t)  :  f(t)  +  A_££t) 


2 

(tT 


e 

The  following  integral  of  this  approximation  is  used  as  an 
approximation  for  M(t) , 

t 
E(t)  -  F(t)  +  X  f  jjl^  du.  (11) 

o 

Another  form  of  the  renewal  function  can  be  derived  from  (1) 

or  (2) .   One  can  show  that 


/ 


[1  -  Fe(t-u)]  dM(u)    Fe(t). 
"o 

If  we  approximate  dM(u)  on  the  right  hand  side  by  Adu  we  get 

another  approximation  for  M.   Thus  define, 

l  n 


t 

6(10  =  At  +  A  f  [1  -  Fe(u)]  du  -  Fe(t).      (13) 
o 

Recall  that  A(t)  and  B(t)  gave  simple  lower  and  upper  bounds 

on  M  for  F  restricted  to  NBUE .   However,  (A(t) -At)-*-l  and 

(B(t)  -  At)->-0  as  t-K) .   From  (3)  for  F  non-lattice  it  is  known 

that  (M(t)-At)+k,  so  define 

H(t)  =  At  -  Fe(t)  +  (1+k)  F(t),  (14) 

the  convex  combination  of  A  and  B  with  the  correct  asymptotic 
behavior. 

6 .   Distribution  Tested 

The  following  distributions  with  specified  param- 
eters are  used  to  compare  the  bounds  and  approximations  for 
M(t). 

a.   Uniform  Distribution  (A=l) 

F(t)  =|  if  0  <  t  <  2 

=1  if  2  <  t. 


Then 


t2 


These  give 


Fe(t)  =  t  -  J-  if  0  <  t  <  2 

=1  if  2  <  t. 


A(t)  =  ^j  if  0  <  t  <  2 

=  t  -  1  if  0  <  t, 

B(t)  =  ^  if  0  <  t  <  2 

=  t  -  1  if  2  <  t, 


11 


coo  ■  i- 


|1+  £ 

4  2 


i£  0    <    t    <   2 


=    t 


if  2   <  t, 


doo  -  ^r 


=  t 


E(t)    =    41og    (  jrt   )    "   I 


=   t    -    41og    (i)-3 


if   0  <  t   <   2 

if  2  <  t, 

if   0  <  t   <    2 

if  2  <  t, 


t3        t2 


G^  =  i2  "  r- +  t 


t-i 


t 2  t 

H(tl   =  J-  +  I 

4  o 


=   *    -    3 


For   this    case, 


M(t)    = 


r  j-i 

Zl       ,      t.n 
EI   ^"^      e 
L  n=0 


!- 


if   0  <    t   <   2 

if  2  <  t, 

if  0  <   t   <    2 

if   2  <   t. 


-1    if    (j-l)(2)    <    t    <    (j)(2) 

™~       X    j   i^    j    •    •    •  • 


b.      Gamma   Distribution    (A=l) 


F(t)    =    1    -    (1    +    2t)    e 


2t 


Fe(t)    =    1    -    (1    +    t)    e"2t. 


These   give 


A(t)    =   t   +    (1    +    t)    e"2t    -    1, 


B(t)    = 


1    -    (1    +    t)    e 


-2t 


1, 


C(t)  =  t  -  t  e'2t, 


D(t)  = 


t  -  t  e 


2t 


2  t2  e 


2t 


1  -  e 


^2t 


t  e 


rTt 


E(t)  =  1  -  e~2t  -  2  t  e"2 


••;^t 


-2u  ~    -2lu  2 
e    -2u  e    ) 

-  2u      -2u 
e    -  u  e 


du, 


G(t)  =  t  +  '(  \   t  +  i  )  e"2t  -  i 
H(t)  =  t  -  |  (  t  -  |  )  e'2t  -  J-. 


For  this  case 


M(t)  =  t  +  \   e"4t  -  J.. 

c.   Hyperexponential  Distribution  (A=l) 

F(t)  =  |  (  1  -  e"2t  )  +  \   (  1  -  e  5  ^  ), 
Fe(t)  =  |  (  1  -  e'2t  )  +  |  (  1  -  e  5    ). 


These  give: 


-2  t 

A(t)  =  t  +  |  e"2t  *  |  e  5    -  1, 


B(t)  = 


8t 


-1  t 
3(1  -  e"2t  )  +  5(1  -  e  5    ) 


C(t)  -  t  +  |   (  1  -  e"2t)  -  |  (  l  -  e  5    ), 


D(t)  = 


2t[3(l  -  e"2t  )  +  (  1 


-ft 

e  5    U 


3(  1  -  e"2t  )  +  5  (  1 


e  5  '  ] 


-It 

E(t)    =|(1-   e"2t   )    +   \   {   1    -   e   5        ) 


Jiia 


_2 
e"2u   )    +    (i    -    e"5   U)]2 


du    , 


u 


0      6(1        e'2u)    +    10(1    -    e    5 


-2   t 
G(t)    =    t    -    -I    (    1    -    e'2t)    +    If    (    1    -    e    5      ), 


16 


16 


.1  t 
H(t)    =   t   +   1|   (   l    -   e"2t)    -   ^   (   1    -    e    5        ) 


For   this    case 


M(t)    =    t   +   |   (   1    -   e    5      ) 


d.   Log-normal  Distribution  (A=/e— ) 
Let 


*(t)  = 


m    J 


t    X: 


e     dx  . 


Then 


F(t)  =  $  (log  t) 


F-Ct)  =  te    (1  -  *(log  t))  +  $  ((log  t)  -  1) 


These  give 


A(t)  =  e  "   t  <D(log  t)  -  *((log  t)  -  1), 

1 


B(t)  =  — T 


e  2  t 


e  L   t(  1  -  »(log  t))  +  $  ((log  t)  -  1) 


1, 


C(t)  =  (e  2  t  +  1)  *(log  t)  -  *((log  t)  -  1), 


D(t)  =  — e  2  t  *ClQg  $ 

e  J   t(l  -*(log  t))  +  $((log  t)  -  1) 


E(t)  =  *(log  t)  +  e 


t 

($(log  u))2  du 


'In 


o    e  zu(l-$(log  u))+$((log  u)-l) 


G(t)  =  e  2  t  -  I   (1  -  *((log  t)  -  2)) 


+  (e  2  t  +  1)  (1  -  •((log  t)  -  1)) 


_1 


-  \   e  2  t  (e  2  t  +  2)(1  -  »(log  t))  +  %   -  1, 


_1 
H(t)  =  (e  2  t  +  I  e)  »(log  t)  -  »((log  t)  -  1) 

For  this  case  the  exact  form  of  M(t)  cannot  be  found. 


III.   NUMERICAL  RESULTS 

Numerical  calculations  are  presented  for  each  bound  and 
approximation  and  all  four  distributions  in  Section  II.   Re- 
sults are  tabulated  below  and  are  graphed  and  comparisons 
are  made  with  M(t) .   For  the  log-normal  case  M(t)  was  simu- 
lated for  comparison  with  approximations. 
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FIGURE  7 


I  2  3  4  5  t 

:  M(t)  and  E(t)  for  the  hyperexponen  al  distribution. 
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FIGURE  8 


:  M(t)  and  G(t)  for  the  hyperexponential  distribution. 


27 


FIGURE 


RE  9  :  M(t)  and  H(t)  for  the  hyper exponential  distribution. 
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